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Abstract. 

Asymptotical formulas for the photon distribution function of a quantum oscillator 
with time-dependent frequency and damping coefficients, interacting with a thermal 
reservoir, are derived in the case of a large mean number of quanta. Different regimes 
of excitation of an initial thermal state with an arbitrary temperature are considered. 
New formulas are used to predict the statistical properties of the electromagnetic field 
created in the experiments on the Dynamical Casimir Effect which are now under 
preparation. 

PACS numbers: 42.50.Ar, 42.50.Lc 
1. Introduction 

The so-called Dynamical Casimir Effect (DCE), i.e., a generation of photons from 
vacuum due to the motion of neutral boundaries, was a subject of numerous theoretical 
studies for almost 40 years (for the specific reviews on the DCE see [IH3]; general 
reviews on the Casimir physics, including different manifestations of the static and 
dynamical Casimir effects, can be found, e.g., in [HE]). One of the most important 
results obtained during a decade before and after 2000 was a prediction of a possibility 
of an observation of the DCE in a laboratory, using high-Q cavities with dimensions of 
the order of few centimeters and resonance frequencies of the order of few GHz or higher. 
Namely, calculations performed by several groups of authors [6T4T2] (using quite different 
approaches) resulted in the same conclusion: if one could arrange periodical changes of 
parameters of some cavity (its dimensions or properties of the walls) for a sufficiently 
long time, then initial vacuum or thermal fluctuations of the electromagnetic field could 
be amplified to a detectable level due to the effect of parametric resonance. This result 
stimulated the work of several experimental groups, from which the MIR group of the 
university of Padua [13] seems to be close to a success (the name of the experiment 
MIR means 'Motion Induced Radiation'; this term was coined in [8]). Therefore it 
seems interesting to calculate the statistical properties of quantum states which could 
be obtained under realistic experimental conditions. 
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The simplest model of the DCE in cavities is that of a quantum oscillator with a 
time-dependent frequency, describing the selected field mode which is in resonance with 
time variations of the cavity parameters. Such a model was first proposed in [r^[T5] 
and later it was developed in [6j[7]. In the ideal case (without losses or interactions 
with other degrees of freedom) the oscillator goes from the initial ground state to the 
vacuum squeezed state [Bj, whose properties are well known. In particular, this state 
has a strongly oscillating photon distribution function (abbreviated hereafter as PDF) 
f(m) = (m\p\m) (i.e., the probability to detect m quanta in the state described by the 
statistical operator p): 

(n) m (2mV 

/(2m) = ,\' V /o ; — , f(2m + l) = 0, (1) 

where (n) is the mean number of quanta. Another important property of the vacuum 
squeezed state is a high degree of quadrature squeezing when (n) 3> 1. 

A quite different situation arises if the intermode interaction is essential (as happens 
in the effectively one-dimensional Fabry-Perot cavity [THE] or i n the three- and two- 
dimensional cavities with accidental degeneracies of the eigenmode spectra [TQlITH] ) or 
coupling with a detector is strong enough |6j. Then the degree of squeezing in each 
mode becomes much smaller than in the vacuum squeezed state and the oscillations of 
the PDF disappear. 

New features appear when the dissipation becomes important. Just such a situation 
takes place in the MIR experiment. Its main idea is to simulate a motion of one of 
the cavity walls using an effective electron-hole 'plasma mirror', created periodically 
on the surface of a semiconductor slab (attached to the wall) by illuminating it 
with a sequence of short laser pulses. If the interval between pulses exceeds the 
recombination time of carriers in the semiconductor, a highly conducting layer will 
periodically appear and disappear on the surface of the semiconductor film, thus 
simulating periodical displacements of the boundary. In this way, rather big relative 
amplitudes of displacements and changes of the fundamental cavity eigenfrequency, of 
the order of 10~ 3 or 10~ 2 , can be easily achieved using standard semiconductor plates 
having the thickness of the order of 1 mm. This is is a great advantage over the schemes 
where real oscillations of the surface of the cavity walls are excited, because in the 
latter case the relative amplitude of displacements cannot exceed the value 10~ 8 due to 
tremendous internal stresses arising inside the material for the frequencies of the order 
of a few GHz or higher. 

Note that the thickness of the photo-excited conducting layer nearby the surface 
of the semiconductor slab is much smaller than the thickness of the slab itself. It is 
determined mainly by the absorption coefficient of the laser radiation, so it is about 
few micrometers or less, depending on the laser wavelength. Therefore laser pulses with 
the surface energy density about few /xJ/cm 2 can create a highly conducting layer with 
the carrier concentration exceeding 10 17 cm -3 , which gives rise to an almost maximal 
possible change of the cavity eigenfrequency for the given geometry [15] . It is worth 
noting that although the thickness of the conducting layer is less than the skin depth, it 
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gives the same frequency shift as the conductor filling in all the slab. This interesting fact 
was explained and verified experimentally in [20] . On the other hand, the conductivity of 
the layer is not extremely high due to a moderate value of the mobility in the available 
materials (such as highly doped GaAs), which is less than lm 2 V _1 s [13]. For this 
reason, effects of dissipation during the excitation-recombination process inside the layer 
cannot be neglected, since they can change the picture drastically [21 [19]. A model 
taking into account the dissipation was developed in [2], and some of its consequences 
with respect to the photon statistics and PDF were considered recently in PET]. Here I 
analyze some interesting special cases that were not considered in the previous papers. 

2. Exact formulas for the PDF and statistical moments 

I assume that the effects of dissipation can be taken into account by means of a model 
based on the Heisenberg-Langevin operator equations of the form (h = 1) 

dx/dt = p- lx (t)x + F x (t), (2) 

dp/dt = -y p (t)p - u 2 {t)x + F p (t), (3) 

where x and p are dimensionless quadrature operators of the selected mode of the EM 
field. These operators are normalized in such a way that the mean number of photons 
equals M = \{p 2 + x 2 — 1). Two noise operators F x {t) and F p {f) with zero mean 
values (commuting with x and p) are necessary to preserve the canonical commutator 
[x (t),p(t)} = i. These operators give a simplified description of complicated processes 
inside a thin lossy dielectric (semiconductor) slab attached to one of the cavity walls. In 
the phenomenological model used in this paper, the net result of all those processes is 
encoded in the correlators of the noise operators (the Markov approximation is assumed) 

(F j {t)F k (t')}=6(t-t') Xj k{t), j,k = x,p. (4) 

In principle, the functions Xjk(t), as well as the functions j x (t) and j p (t), should 
be derived from some 'microscopical' model of electron-photon and electron-phonon 
interactions inside the semiconductor slab. In the phenomenological model considered 
here, these coefficients can only be 'guessed'. The simplest possibility consistent with 
principles of quantum mechanics is as follows [21I2T]: 

lx (t) = 7p (t)= 7 (t), (5) 
Xx P = -X P x = ij(t), Xxx = X PP = l(t)G, (6) 

where 7 (i) is identified with the imaginary part of the complex time-dependent 
eigenfrequency of the cavity uj c (t) = oj(t) — z 7 (t), which can be found from the solution 
of the classical electrodynamical problem by taking the instantaneous geometry and 
material properties (for example, by solving the Helmholtz equation with the complex 
dielectric function e(t) inside the slab and the set of cavity dimensions {Lj(t)}, where 
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the time variable t is considered as a parameter). The coefficient G is related to the 
temperature of the reservoir as 

G = 1 + 2(n) th = coth [hwi/(2k B e)] (7) 

(so that (n) t h is the mean number of quanta for the given cavity mode in the 
thermodynamic equilibrium). The specific set of coefficients (jSJ) and ((6]) has the following 
remarkable property: if the frequency u does not depend on time, then the second-order 
statistical moments (x 2 ), (p 2 ) and (xp + px) tend to the well known equilibrium values 
for an arbitrary positive function 7(£) (without any corrections). Some arguments pro 
e contra the choice of damping coefficients in the form (j5]) can be found in |2T]. 

Equations (J2J) and (EJ) can be solved explicitly for arbitrary time- dependent 
functions ■y x , P (t), co(t) and F x>p (t): 

x(t) = x s (t)+X(t), pit) = p s (t) + P(t). (8) 

The first terms are the solutions of homogeneous equations (without the noise operators) 

x s (t) = e~ r « {x Re [£(£)] - p lm [£(£)]} , (9) 



Psif) = e~ r « {x Re [i(t)] - polm [£(t)] } , 



(10) 



where xq and po are the values of operators at t = (taken as the initial instant) and 
r(t) = J*7(r)dT. Function £(t) is a special solution to the classical oscillator equation 

e + u; 2 (t)£ = 0, (11) 

selected by the initial condition £(£) = exp(— it) for £ — > — oo, which is equivalent to 
fixing the value of the Wronskian 

The operators X(t) and P(£) represent the influence of the stochastic forces: 

( x(t) \ 
V Ht) J ' 



-r(t) 



dre T{T) A{t; r) 



(12) 



(13) 



where the 2x2 matrix 
A(t; t) = 



consists of the following elements: 



Im 
Im 



mt(r) 
mt(r) 



<(t;r) \ 

aP v {t-r) J 
al = Im [C{mr)\ 



Im 



(14) 

(15) 
(16) 



Combining equations f|T3l - ffT6l) with (SD-([H]) one can obtain the following expressions 
for the second-order moments of operators P(t) and X(t) [here ft = /(£)]: 



(X 2 (t)} 



\U 2 J t -Re(^ 2 J t ), 
|£ t | 2 J t - Re (ef J t ) , 



(17) 
(18) 
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1 -{XP + PX) t = Re - CtCtJt) , (19) 

Jt = ^- 2rW [dre^h{r) (%\ 2 + |^| 2 ) , (20) 

4 = f e" 2r(t) pre 2r « 7 (r) (£ + g) . (21) 

The mean number of quanta can be written as M{t) = Af s (t) + Af r (t), where the 
first term depends on the initial state ('signal'), while the second term is determined by 
the interaction with the reservoir. From (TPTj) and (|18p one obtains 

K(t) = E t J t -Re(E*J t ), (22) 

E t =l(\tt\ 2 + \ti 2 ), Sk = l(& + H)- (23) 

For the initial thermal state characterized by the parameter Go (which can be different 
from G) one has 

Af( th \t) = l - {G e- 2T ^E(t) - l} . (24) 

One should remember that formulas (122!) and (1241) make sense for sufficiently big values 
of time t, when the laser pulses have been switched off and the recombination processes 
have been over, so that the normalized frequency u(t) returns to its initial unit value 
(because the photon number operator N — | (p 2 + x 2 — 1) is defined with respect to 
the initial geometry of the cavity, coinciding with the final one). 

I consider here only the special (although the most realistic) case of initial thermal 
states of the field. It is well known [221123] that the description of open quantum systems 
by means of the Heisenberg-Langevin equations with delta-correlated stochastic force 
operators is equivalent to the description in the Schrodinger picture by means of the 
master equation for the statistical operator. In the case of linear operator equations of 
motion, such as equations fl2]) and (j3J), the corresponding master equations contain only 
quadratic terms (various products of two operators p and x) [2"2"H2"5] ). Consequently, 
any initial Gaussian state (whose special case is the thermal state) remains Gaussian 
in the process of evolution. For thermal states, mean values of the first order moments 
are equal to zero, and this property is preserved in the process of evolution governed by 
equations ^ and Then all statistical properties of the single mode are determined 
completely by the variances of the quadrature operators a xx , a vv and by their covariance 
o~x P = o~ px (in the case involved, a a b = h(ab + ba)). Using equations (jUJ) and ( ITU]) one can 
verify that the time- dependent (co) variances can be obtained from formulas f|T7]) - (|T9]) 
by means of the replacement J t — > J t + |G exp(— 2Y t ). 

The photon distribution function f(m) = (m\p\m) of the Gaussian states was found 
long ago [26H3T] . For zero mean values (x) — (p) — it can be expressed in terms of 
the Legendre polynomials 

2D m/2 ( 4A - 1 \ 
ff m ) = ~ /0 P m (25) 

M ) D (m+l)/2 m \^D~fDZ ^ ' 
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where 

D± = 1 + 4A ± 2r, (26) 

t = a xx + a pp = 1 + 2J\f, (27) 

A = <T xx a pp - a 2 px > 1/4 (28) 

(the last inequality is the Schrodinger-Robertson uncertainty relation). For the 
Gaussian quantum states one can write A = l/(4/z 2 ), where the quantity \x is the 
quantum purity of the state: /i = Tr(p 2 ). For the pure Gaussian quantum states (in the 
absence of dissipation) A = 1/4, and formula ( 1231) goes directly to (JTJ) - 

The explicit expression for the time- dependent coefficient A is as follows: 



2 



A =^J+^e- 2r j -\J\ 2 . (29) 

If the functions u(t) = uo[l + xif)} and ^(t) have the form of periodical pulses 
separated by intervals of time with u = Uq = const and 7 = (this means that the 
quality factor of the cavity is supposed to be high enough), then the following formulas 
can be obtained for the quantities E n = E(nt), E n = E(nt), J n = J{nt) and J n = J{nt) 
after n periods (n pulses of laser irradiation in the case of DCE) under the realistic 
conditions \x(t)\ < 1 and >y(t) < 1 0121]: 

E n = cosh(2nz/), E n = sinh(2ni/)e i/3 , (30) 
J n = + \ J n = e i/3 (^A^ — A^j , (31) 

where /3 is some insignificant constant phase, 



v 



co x{t)e- 2iW0t dt 



A= / 7(r)dr. (33) 



Here U and t/ are the initial and final moments of each pulse. It is taken into account 
that A, v 1. The mean number of photons grows exponentially under the conditions 
2nv > 1 and v > A: 

M n = -e 2n ^ (g + -^-) + 0(1). (34) 



4 V u v-k, 

Coefficient Go is given by formula (jTJ), but with 6 replaced by the initial temperature 
of the field mode Oo (which can be made different from B). Under the same conditions 
formula (129]) takes the form 

C A 

A n =A/;— ^ + 0(1)- (35) 

Note that the ratio A n /Af n in this limit does not depend on the coefficient Go (the 
initial temperature of the field mode). Formulas (130|) - (135|) hold provided the periodicity 
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of pulses T is adjusted to the initial period of oscillations of the field mode T = 2n/u 
as follows (m = 1,2,...): 

T = ]-T (m - <p/ir) , <p = -u> [ tf x{t)dt. (36) 

A small shift (p of the resonance periodicity of pulses arises if the profile of pulses is 
asymmetrical (namely this situation takes place in reality). 



3. Asymptotical formulas for the photon statistics and squeezing 

Formula (125]) is exact. However, it is not very convenient for calculations in the case 
of DCE if the number of created photons is big (say, m ~ M > 1000; otherwise the 
effect cannot be confirmed with certainty at the existing experimental level due to the 
noise in the measurement channel). Therefore asymptotical forms of exact formulas for 
m ^> 1 can be more useful. Note that the argument of the Legendre polynomial in ( 1251) 
is always outside the interval (—1, 1), being equal to unity only for thermal states with 
r = 2yA. For this reason it is convenient to use the asymptotical formula [32] 

P m (coshO« (-^) ' /o([m + l/2]0, (37) 

where Iq(z) is the modified Bessel function. Formula (T25]) shows that the behavior of 
the PDF depends on the sign of the coefficient _D_, i.e., on the ratio 2A/r rj A n /J\f n . 
If this ratio exceeds the unit value, then the argument of the Legendre polynomial 
is real (and bigger than unity). In this case, using the known asymptotical formula 
Io(x) ~ (27tx)~ 1//2 exp(x) (if x ^> 1) and making some further simplifications, one can 
arrive (under the conditions J\f ^> 1 and m ^> 1) at the simple formula [3"|12Tj 

w exp[-( m+ l/2)/(2AQ] 

^27rA^(m+ 1/2) 

Here I consider in detail the case when D- < or A n /J\f n < 1. Then formula (13Tj) is still 
valid, but it needs some transformations because the parameter £ becomes complex. It is 
clear from formula (125]) that the real positive value f(m) does not depend on the choice 



of sign of the square root function yTTZ = ±iy |Z?_|, provided this sign is maintaned, 
both inside the argument of the Legendre polynomial and in the coefficient in front of 



this polynomial. Choosing for definiteness the branch yUZ = zy \D-\ it is convenient 
to write 

4A - 1 

£ = £ — 27r/2 = — irj, sinh(£) = 



DADJ 



so that £ is real positive number. Using the relation I^—irf) = Jo{rj) in (137j) . one can 
represent (125*1) as 



f{m) 





f 




V r / 







J Q ([m+ 1/2] 77) 
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where r = \Jt 2 — 4A. Since arg(r^) < ir/2, the asymptotical formula Jo(z) 
— cos(z — 7r/4) can be used, leading to the formula 

2z m cos (mir/2 + i [m + 1/2] f 



/(m) 



wV(m + l/2)r 



Replacing the cosine function by the sum of two exponentials, one obtains 



f{m) 



1 



\Jn(m + l/2)r 

+ (-ir ' ID " 



'2r + <5' 



m+7 



2r + 5, 



(39) 



where 5 = 4 A — 1. The right-hand side of ( 139]) is real and positive. Using the relation 
2 A = br + c for r 3> 1 [with b = GA/ {y + A) according to equation (135|) ]. one can verify 
the relations 

2r + 5 

= 1 - 

T 



l-D- 



(40) 



1 



2r + S 1 + 6 



2c 



r(l -6 2 ) 







(41) 



Note that the coefficient at r" 1 in ( 140]) does not contain the coefficients b and c. 
Obviously, the expansion (jUJ) is valid provided 1 — b is not too small. Formulas ( 1391 ) 
and (jHJ) clearly show that oscillations of the PDF are actually negligible (especially for 
m ~ TV), unless b 1 (i.e., in the case of extremely small dissipation). In the latter 
case the coefficient c can be expressed in terms of the initial purity which, in turn, is 
related to the factor G as ji = 1/G , so that c = Gq/2. Using the approximate formula 
(1 — x) m w exp(— mi), which holds for i C 1 and mx 2 1, one can rewrite (139]) in the 
case of small dissipation as (replacing r by r) 

[7rr(m + 1/2)] ~ 1/2 {exp [-(m + 1/2) /r] 



/(m) 



+ (-l) m exp 



— m 



l/2)G 2 /r]}. 



(42) 



Formula (j4"2l) is valid under the conditions r w 2J\f > 1, 1 Cm Cr 2 and Gkjv <C 1. 
Even in this case the oscillations of the PDF can be noticed only if Go ~ 1, i.e., for low 
temperature initial thermal states. In particular, (142!) coincides with the asymptotical 
form of the ideal PDF of the vacuum squeezed state (OQ) if Go = 1 and r ~ 2(n) 1. In 
the conditions of the MIR experiment the ratio A/u exceeds 1/2 |3], which means that 
the oscillating term in (139]) can be neglected even for G = 1 (zero temperature of the 
cavity walls) , so that formula ( 138]) can be used for any values of parameters Go and G 
(under the conditions Af ^> 1 and m 3> 1). 

Using the Euler-MacLaurin summation formula, one can verify that the distribution 
function (138]) has the correct normalization with an accuracy (9(r -1 / 2 ): 



m=0 



f(m)dm + O[f(0)\ 



(43) 
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eM-xjr)_ dx + 0,-1/2) = l + 0(r -i/2). (44) 

y'KTX 

The oscillating terms in ( 139]) or ()42p do not influence the normalization, because they give 
corrections of a higher order. For example, combining the nearest positive and negative 
terms with m = 2k and m = 2k + 1 and applying the Euler-MacLaurin summation 
formula to the sum over these positive pairs (which contains now only slowly varying 
positive terms), one can see that the correction equals approximately G /(2t). 

With the same accuracy as in ( 1441) . the moments of the distribution function can 
be calculated as 

(m k ) = Ym k f(m)n f°° x fe exp( Zf/ r) dx (45) 

= T k {2k ~ m *Af k (2k-l)\\. (46) 

For k = 2 formula f l4^|) leads to the following formula for the variance of the number 
of created quanta: = (m 2 ) — (m) 2 m 2M 2 . It means that the field mode goes 
asymptotically to the so called 'superchaotic' [33l|3l] quantum state, whose statistics 
is essentially different from the statistics of the initial thermal state, characterized by 
formula = M{M + 1) ~ M 2 . The same result can be obtained from the general 
formula for the variance of the number of quanta in the Gaussian quantum states [30] 
^ = r 2 /2 - A - 1/4 (if ( x } = (p} = 0). 

It is interesting that the asymptotical value of the ratio 2A/r = b in the case 
concerned coincides with the invariant squeezing coefficient. This coefficient can be 
introduced in the following way. Obviously, instantaneous values of variances a xx , a pp 
and a xp cannot serve as true measures of squeezing, since they depend on time in the 
course of the free evolution of the oscillator. For example (in dimensionless units with 
wo = 1), 

a xx (t) = cos 2 (t) + 4°J sm 2 (t) + a x °J sin(2t), (47) 

and it can happen that both variances a xx and large, but nonetheless the 

state is highly squeezed due to the large nonzero covariance a xp . It is reasonable to 
introduce some invariant characteristics which do not depend on time in the course 
of free evolution (or on phase angle in the definition of the field quadrature as 
E((p) = aexp(— iip) + a) exp(iy?) /a/2). I define the invariant squeezing coefficient 
5* as the ratio of the minimal value of the variance a xx (t) [as a function of time ( 14T|) ] to 
the dimensionless variance 1/2 in the vacuum state. Then straightforward calculations 
give the formula (similar results can be found in [3TII351437] ) 

S = 4A (48) 

t + Vt 2 -AA 

so that, indeed, S = b = GA/(u + A) in the case discussed. In contradistinction to the 
ideal case (7 = 0), when A = const, so that an arbitrarily small invariant squeezing 
coefficient Sid ~ G\l (4A/") can be obtained, in the presence of dissipation the squeezing 
coefficient goes to the constant value which does not depend on the initial temperature 
of the mode but depends on the temperature of the reservoir (the cavity walls). 
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4. Conclusion 

The main results of the paper are as follows. New formulas ( 1391) and (142 1) show how 
the increase of dissipation and temperature transform the strongly oscillating photon 
distribution function (TI]) of the ideal squeezed vacuum state to the smooth asymptotical 
distribution ( 1381) . Namely this smooth distribution is expected to be observed in the MIR 
experiment due to strong dissipation. Another new result is the explicit demonstration 
of the correlation between the existence of squeezing (when the invariant squeezing 
coefficient S is less than unity) and the oscillations of the PDF in the case concerned. 
In particular, no squeezing of the fundamental field mode is expected in the MIR 
experiment. 
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